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Abstract 


Buoyancy-driven flow in a superimposed air-water layer has been experimen- 
tally studied, using a Mach-Zehnder interferometer. Two methodologies were 
employed, namely experimental and numerical. The cavity containing the fluid 
layers has dimensions of 447 x 32.1 x 32.1 mm^. The cavity temperature differ- 
ence was maintained at AT=10°C during the experiments. The fluid properties 
and Rayleigh numbers have been evaluated at the average cavity temperature. 
The sequence of experiments conducted were with a (i) cavity filled with air, 
(ii) cavity filled with water and (iii) cavity half-filled with water. All fringe pat- 
terns have been recorded in the infinite fringe setting of the interferometer. The 
experiments have been conducted till steady state is reached. 

To examine the experimental results, a two-dimensional numerical simulation 
has been carried out for flow and heat transfer in air (Pr=0.71) and liquid metal 
(Pr=0.05). Closed as well as open cavities have been considered. 

The present research is directed towards examining the following issues: (a) 
the influence of Rayleigh number on the steady thermal field, (b) the transient 
evolution of the thermal field at a given Rayleigh number, (c) temperature dis- 
tribution at the interface, (d) heat flux at the solid wall, and (e) the shape of the 
interface 

In air, a steady state two-dimensional convection pattern was observed. The 
average Nusselt number for the cavity matched well with the published correla- 
tions. In water filled cavity, due to an extremely high Rayleigh number, unsteady 
and turbulent flow was observed. For a cavity half-filled with water, the inter- 
face temperature calculated from the water side was found to be quite close to 
that from the air side. The direction of circulation in the air phase was found 
to depend on the average cavity temperature. For a cavity average temperature 
of 21°C, mechanical coupling was observed between the two layers and for 29°C, 
thermal coupling was observed. 
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Chapter 1 
Introduction 


The heat transport processes augmented by fluid motion in a horizontal fluid 
layer confined between a pair of parallel plates heated from below and cooled from 
the above is often called Rayleigh-Benard Convection. Due to the temperature 
gradient imposed across the plates, the fluid layer near the lower wall expands and 
thereby becomes less dense than the overlying layers. This is a gravitationally 
unstable configuration. In order to achieve stability, the warmer lighter layer at 
the bottom will tend to rise and the cooler heavier layer at the top will tend to 
sink. This will result in a convective flow, leading to a cellular pattern of fluid 
motion. 


1.1 Convection in Fluids 

In Rayleigh-Benard convection a thin layer of the fluid is confined between two 
flat, rigid, horizontal plates. The fluid completely fills the space between the 
plates. In particular, there is no free surface formed. The fluid layer is heated 
from below and cooled from the top resulting in a convective motion. Buoyancy 
is the driving force responsible for the convective motion. If buoyancy force is 
smaller than a certain threshold, the viscous forces between the fluid layers acts 
as internal friction and inhibit fluid motion. Heat transfer across the fluid layer 
is then purely by conduction. Once buoyancy forces exceed the threshold value, 
convective motion starts. 

The non-dimensional parameter Rayleigh number Ra is linked with buoyancy 
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and viscous forces in terms of energy as 

Ra — potential energy _ gP{Thot - Tcoid)H^ 

viscous dissipation ua 

The Prandtl number is a second dimensionless quantity arising in thermal con- 
vection. It is a measure of ratio of the molecular diffusivity of momentum to that 
of thermal energy and is a fluid property. The Prandtl number is defined as 


p _ molecular diffusivity of momentum u , . 

molecular diffusivity of thermal energy a 

Flow transitions are generally described in terms of Rayleigh number and Prandtl 
number. 

Two-layer convection studied in 'the present work comprises of differentially 
heated superposed fluid layers. The layers are confined between two solid, hori- 
zontal isothermal surfaces heated from below and cooled from top. It is a trans- 
port problem that is affected by the ratios of thermo-physical properties such 
as ji, p, K of the two fluid layers and geometrical factors. The main difference 
between one and two-layer convection is in the formation of a free surface. This is 
possible with a liquid-gas combination which gives rise to surface-tension-driven 
convection apart from being buoyancy-driven. In superimposed immiscible liq- 
uids, two-layer convection gives rise to deformable interfaces. Interfacial surface- 
tension-gradient-driven convection is caused by the variation of surface tension 
with temperature. Like density, the surface tension decreases as temperature 
increases. Therefore the cooler region will tend to draw fluid particles towards 
itself in comparison to warmer region. 

Marangoni number Ma is a non-dimensional parameter which links surface 
tension and viscous forces in terms of energy as 

,, surface tension energy cttATH 

Ma = — ^ = (1.3) 

VISCOUS dissipation /tk 


Convective motion can be initiated by surface tension alone, even in the absence 
of a gravity field. 

Benard-Marangoni convection combines the effect of buoyancy and surface 
tension. This typically occurs in a liquid-gas system. Buoyancy is initiated in 
the upper gas layer, and simultaneously sends a thermal signature to the surface. 
This generates either surface tension or buoyancy convection in the lower layer. 
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Rayleigh-Benard-Marangoni convection is of interest in a variety of engineer- 
ing applications and has a great scientific importance. It has direct relevance to 
atmospheric convection. The importance can be seen from geophysical studies to 
astronomical research. The need to study fluid convection is very strongly felt 
for theoretical as well as fundamental reasons. 

The Rayleigh-Benard convection is an important mechanism of fluid flow and 
heat transfer in natural as well as many engineering systems. Important applica- 
tions of two-layer convection are: 

1. convection in the earth’s mantle. 

2. liquid-encapsulated crystal growth. 

3. development of new generation boilers. 

4. weather prediction programme. 

5. glass processing in low gravity conditions. 

1.2 Literature Review based on the Two Layer 
Convection 

Literature for two layer convection can be categorized into: 

1. Analytical studies 

2. Numerical studies 

3. Experimental studies 

Convection in these studies could be dominated by buoyancy alone, surface 
tension alone or a combination of both. 

1.2.1 Analytical studies 

Zeren and Reynolds (1972) predicted that the instability response is dependent 
on the properties of the fluids. Instability also varies with the total depth of layer 
and depth fraction of one fluid. 
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Through a mathematical model Busse (1989) has investigated the onset of 
convection in two layers. The two superimposed immisible fluids with combina- 
tions of ethylene glycol-oil and ethylene glycol-decane were heated from below. 
He showed that the convection in two layers may occur in the form of either 
viscous or thermally coupled motions. He concluded that oscillatory instability 
exists and corresponds to a cyclic variation between above motions. Narayanan 
(1996) has determined the pattern formation at the bifurcation point from the 
quiescent state. The critical conditions which determine the critical Marangoni 
number Ma for different aspect ratios were found. The patterns for the onset of 
three-dimensional convection in a bounded circular container were theoretically 
obtained. Golovin (1997) presented a theory of secondary instabilities in silicone- 
oil system with a deformable interface. He concluded that this instability is seen 
to be either monotonic or oscillatory. 

1.2.2 Numerical studies 

Busse (1996) has analyzed the complex interaction between rolls of different ori- 
entations. The stability of undulating roll states and their modulations in the 
two-layer system were also analyzed. He concluded that through controlled initial 
conditions rolls with wavelength in the zigzag unstable regime can be generated. 

Johnson and Narayanan (1997) concluded that if depths of fluids layers are 
approximately equal, then the respective convective fields can have two couplings: 
viscous (or mechanical) or thermally. In viscous coupling, the convection rolls 
rotate in the opposite direction while in thermal coupling, they rotate in the 
same direction. Dijkstra (1992) studied the problem of pattern selection in small 
aspect ratio containers. Because of limited computational resources, he restricted 
the work to two-dimensional containers. He concluded that no-slip sidewalls 
greatly influence the multiplicity of the stable steady patterns. Nield (1964) has 
combined the two rival theories of buoyancy-driven and surface tension-driven 
convection. He showed that the two theories causing instability reinforce each 
other and are tightly coupled. The corresponding cells formed are of same size. 
Dauby et al. (1997) presented a study of Benard-Marangoni instability in small 
rigid container. They determined the critical conditions for onset of motion as 
well as the convection patterns at the threshold. 
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1.2.3 Experimental studies 

Andereck et al. (1996) have presented preliminary results from their experiments 
in a system consisting of two immiscible fluid layers in a rectangular geometry, 
driven by a vertical temperature gradient. Time-dependent variations in the 
type of coupling observed between layers have been characterized and described. 
Prakash and Koster (1996) have studied the two-dimensional convection in a two 
layer system the with real-time holographic interferometry. They concluded that 
thermal and mechanical coupling mechanisms are both possible between the lay- 
ers. Johnson and Narayanan (1996) studied the appearance of the codimension- 
two point in a siliconeoil-air combination placed between conductive plates. The 
upper plate consisted of zinc selenide, which is transparent to infrared radiation. 
The flow visualization technique used an infra-red IR camera as the sensor. It 
was placed above the apparatus and viewed down upon the surface to catch the 
flow patterns at the interface. 

Schatz et al. (1999) recently investigated the transition between hexagonal 
and square patterns. It occured as the heating is increased in a fluid of Prandtl 
number 81. The flow visualization technique used in the silicone-air system was 
Shadowgraph. Busse (1996) has discussed two-layer convection from experimental 
as well as theoretical point of view. He presented the results from an experiment 
with two-fluid layer in an annular channel with 1:2 ratio of layer heights. 

Johnson, Narayanan, and Dauby(1999) in their experimental studies have 
presented the effect of air height on the pattern formation in liquid-air bilayer 
convection. For visualizing the flow patterns in the silicone oil-air combination 
in the test cell an inframetric camera was used. The camera was directly placed 
above the test- section. The camera measured the infrared radiation being emit- 
ted by the silicone oil and thus helped in finding the temperature distribution at 
the interface. 

1.3 Objectives of the present study 

The present study is concerned with Rayleigh-Benard convection in superposed 
fluid layers of air and water. The cavity is rectangular in plan and square in 
cross-section. The length of the cavity is 447 mm and the height is 32.1 mm. 
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Five different experiments are done. In these experiments, the cavity tempera- 
ture difference was kept constant at AT=10°C. The values of Rayleigh numbers 
have been calculated by using Equation 1.1 and the fluid properties have been 
evaluated at the average cavity temperature in all the experiments conducted. 

The present work is directed towards achieving the following objectives: 

1. The influence of Rayleigh number on the steady state thermal field. 

2. The transient evolution of the thermal field at a given Rayleigh number. 

3. Temperature distribution at the interface. 

4. Heat flux at the solid walls. 

5. Shape of the interface. 

It is proposed to study buoyancy-driven convection in superposed layers of air 
and water using laser-interferometry. The fluid layers are taken to be enclosed in 
a cavity. The interferograms have been obtained as a collection of fringe patterns. 
The interfrograms give line-averaged information along the direction of the light 
ray. The interferograms can be used to identify the flow pattern prevailing within 
the cavity. 



Chapter 2 

Apparatus and Instrumentation 


A Rayleigh-Benard experimental setup is apparently simple in design. It com- 
prises of two horizontal surfaces of high thermal conductivity which are main- 
tained at different constant temperatures. The vertical side walls are perfect 
insulators. Though simple in design, various factors are taken into account dur- 
ing the fabrication of the test cell and also while performing the experiments. 
For example, uniformity and constancy of surface temperatures, parallelism of 
the walls defining the fluid layer, parallelism of the optical windows and prop- 
erties of the insulating surfaces are all factors that determine the quality of the 
experiments. Further, some extraneous factors such as building vibrations, air 
motion and changes in the ambient temperature strongly affect the recorded in- 
terferograms. Hence these experiments have to be performed with due care and 
precautions. 

The apparatus used to study the convection phenomenon in the horizontal 
layer of fluid is shown schematically in Figure 2.1. The test cell consists of three 
sections namely the top tank, the test section and the bottom tank. The effective 
length of the fluid layer is 447 mm. The fluid layer is confined by two aluminium 
plates of thickness 3 mm above and below. The two longer sides are confined by 
insulated walls and the other two square sides by optical windows. The hot and 
the cold surfaces were maintained at uniform temperatures by circulating water 
over them. To reduce edge effects, the heat transfer areas of the two tanks were 
made quite large compared to the cavity. 

For the designing of the test cell for Rayleigh-Benard convection the important 
steps taken are as follows: Firstly air and water were selected. 



Constant temperature 
I bath (cold) ^ ^ 
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Figure 2.1: Schematic of the Test Cell 
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The selection was based on the transparency of two fluids for studying bi- 
layer convection using laser interferometry. The thermophysical properties of 
fluids were selected in a specific temperature range that meets the experimental 
requirement. 

The law of conservation of energy holds here while dealing with heat transport 
in fluids. Thus heat conducted to water from the lower wall is to be equal to heat 
convected to air which in turn is conducted to the top wall. With this approach, 
the interface temperature is calculated using the Nusselt number correlation given 
by Hollands et al. (1975). The correlation for air is 

N. = l + x.44[l-iMl + l(||)V3-l, (2.1) 


For water, the correlation is given by: 

Nu = 1 + 1,44[1 - 11^1 + (( _ 1 , + 2(^[1 - („<^1 ( 2 . 2 ) 

After getting the interface temperature Tj, and layer heights H, the Rayleigh 
number for each of the fluids is evaluated for air and water by using Equation 1.1. 


Experimentally recorded regimes in Rayleigh-Benard convection (after R.Krishnamurti, 1970) 


TURBULENT 
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3 


X 
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STEADY 2-DlMENSIONAL FLOW 
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100 
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1000 


10000 


Figure 2.2: Various flow regimes in a horizontal fluid layer heated from below 
The dimensions of the cavity which is square in cross-section, have been fl- 
nalised on the basis of individual Rayleigh numbers. The length of the test 
cell was decided on the basis of aspect ratio A. To realize an infinitely long two- 
dimensional cavity aspect ratio should be high. The height of the cavity is decided 
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on the basis of beam diameter (70mm), size of the optical window (50mm) and 
the size of the rectangular tube commercially available. 

The individual heights of air and water have been selected on the basis of 
expected nature of convective flows eis per Krishnamurti’s (1970) flow regime 
diagram, Figure 2.2. In the present study the respective air and water heights 
were used to get the steady two-dimensional convective flow in each layer. 

After getting the dimensions of the cavity, the material is selected for the 
top and bottom walls. Aluminium owing to its high thermal conductivity was 
selected for the two walls. Aluminium provides isothermal conditions at the top 
and bottom walls. Thus perfectly conducting walls required for Rayleigh-Benard 
convection can be achieved. The thickness of the wall were selected by considering 
the two factors (a) minimum thickness for negligible temperature gradient and 
(b) sufficiently thick walls to support water tanks. Thus optimally 3 mm thick 
wall was selected for forming the cavity. 

Now the spacer was chosen on the basis of two factors (a) providing accurate 
vertical spacing between two walls and (b) acts as an insulated side walls. Highly 
polished polyglass (perspex) with very low thermal conductivity was chosen to 
meet the above requirements. 

For the visualization of flow the high grade optical windows of material Fused 
Silica of thickness 10 mm was selected for the front and back walls of the test 
cell. The term high grade refers to the excellent surface quality and parallelism 
of optical window surfaces. 

2.1 Laser Interferometry 

For the temperature measurement in the fluid medium, the primary instrument 
employed in the present work is the Mach-Zehnder interferometer. Figure 2.3 is 
a schematic diagram of the Mach-Zehnder interferometer. 

The optical components of the interferometer are the beam splitters BSl and 
BS2, mirrors Ml and M2. These components are oriented at 45° angle with 
respect to the laser beam direction. The first beam-splitter BSl splits the incoming 
collimated beam into two equal parts, the transmitted and the reflected beams. 
The transmitted beam is the test beam whereas the reflected one is the reference 
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beam. The test beam passes through the experimental apparatus where the 
convection process is in progress. This test beam is reflected by mirror Ml and 
recombines with the reference beam on the plane of second beam splitter BS2. 
The reference beam undergoes a reflection at M2 and passes through the reference 
medium and is superimposed with the test beam at BS2. 

Due to superposition of the test and the reference beams, an interference pat- 
tern is produced at BS2. This pattern contains the information of the variation of 
the refractive index in the test region. For measurements in liquid a compensation 
chamber is required to introduce an appropriate reference environment. 

The mirrors and beam splitters used in this conflguration are of 150 mm 
diameter. Beam splitter has 50% reflectivity and 50% transmitivity. The mirrors 
are coated with 99.9% pure silver and employ a silicon dioxide layer as a protective 
layer against oxidation. 

The Mach-Zehnder interferometer can be operated in two modes; namely 
(a) Infinite fringe setting and (b) Wedge fringe setting. In (a) the test and 
reference beams are set to have zero path difference in the absence of any thermal 
disturbance and fringes form due to density and temperature changes alone. Since 
each fringe is a line of constant phase, it is also a line of constant refractive 
index, a line of constant density and temperature and hence an isotherm. The 
infinite fringe setting is employed for high accuracy temperature measurements 
in the fluid. In (b), the mirrors and beam splitters are deliberately misaligned to 
produce an initial fringe pattern of straight lines. When a thermal disturbance 
is introduced in the path of the test beam, these lines deform and represent 
temperature profiles in the fluid. In general for the measurement of wall heat 
flux the wedge fringe setting is preferred. 

A 35 mW, continuous wave He-Ne laser (A=632.8 nm) is employed as the 
coherent light source for the interferometer. The original laser beam is of the 
diameter 2 mm. A spatial filter is used here to expand the beam to 70 mm 
diameter and a planoconvex lens is used to collimate the beam. 

A CCD (Charge Coupled Device) camera of spatial resolution of 512x512 
pixels has been used to capture the interferograms. This camera is connected 
to a PC-based image processing system through an 8-bit A/D card. The fringe 
pattern is stored in the form of a matrix of integers, with intensities varying 
between 0 and 255 , where 0 indicates black and white is 255. In the frame 
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grabber used, the image acquisition is at the rate of 50 images per second. 

The thermocouples which are used to monitor the temperature of the two 
horizontal aluminium surfaces of the cavity and the ambient temperature were 
connected to a 30 channel recorder. The temperature of the plates were monitored 
throughout the experiment. In the experiments 5 K-type thermocouples are used. 

The optical components of the interferometer are extremely sensitive to vi- 
brations. To avoid ground vibration from reaching the optics, the entire inter- 
ferometer is placed over four pneumatic isolation mounts. These mounts are 
connected to an air compressor for pressurization. Once the mounts are pressur- 
ized the entire interferometer floats over the mounts. The operating pressure for 
the mounts is 5 kg/cm^. This stabilizes the interferometric images and facilitates 
image acquisition. 

A traversing mechanism is needed to mount the apparatus in which the de- 
sired experiment is in progress. This enables translation and rotation of the 
test cell with respect to the flxed positions of the laser source and the optics. 
The traversing mechanism is schematically shown in Figure 2.4. The base of the 
traversing mechanism is padded with a rubber sheet of 30 mm thickness to damp 
any external vibration from reaching the test cell. 

2.2 Alignment of the Interferometer 

Before the start of the experiment the interferometer has to be aligned. All the 
experiments in the present work are performed with the inflnite fringe setting 
mode of the interferometer. Adjustment of the infinite fringe setting is delicate 
and requires efforts. The initial field of view is one of complete brightness since the 
interference is constructive. The geometrical and the optical path lengths of the 
test and reference beams are same in the absence of any thermal disturbance in the 
path of the test beam. As the interferometer approaches the infinite fringe setting 
the distance between the fringes increases and the number of fringes decreases. 
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Figure 2.5 shows the infinite fringe setting mode of the interferometer. 



Figure 2.5; Infinite Fringe Setting of the Interferometer 

2.3 Experimental Procedure 

Firstly the pneumatic isolation mounts, on which the interferometer lies are pres- 
surized. This ensures that the fioor vibrations and other external disturbances 
do not disturb the initial setting. Before the start of the experiments, the Mach- 
Zehnder interferometer is set in the infinite fringe setting mode. 

Isothermal conditions on the lower and upper plates of the test cell are ensured 
by circulating a large volume of water in the lower and upper tank. Temperature 
measurements are made at regular intervals of time during the experiments. Total 
experimental run time of around one and half hours enables the plates to acquire 
a constant and uniform temperature through out their surfaces and hence leading 
to isothermal conditions. 

The interferograms have been recorded at regular intervals of time in order 
to show the transients that the fluid undergoes before reaching the steady state. 
At steady state, three to four interferograms have been recorded in a gap of 10 
minutes to ensure that the steady state condition prevails in the cavity. In all 
the experiments the position of the CCD camera is fixed. The lights of the room 
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are switched-ofF at the time of recording of the interferograms. 

The complete view of the experimental setup with the accessories are shown 
in Figure 2.6 and figure 2.7 shows the closer view of the test cell used to study 
the convection patterns of two phase fluid layer. 



Figure 2.6: Complete View of the Experimental Setup 



Figure 2.7: Closer View of the Test Cell 
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2.4 Uncertainty and Measurement errors 

Errors in the experimental data are associated with misalignment of the apparatus 
with respect to the light beam, image processing operations including filtering, 
thinning and assigning temperature of fringes. Errors related to refraction effects 
in case of water is found to be high. All experiments were conducted several times 
to establish the repeatability of the fringe patterns. In the two phase problem, 
main error occurs during the experiment in air phase due to the dew formation 
on the optical windows. For ail experiments except the one in which cavity is 
fully filled with water, the plate-averaged Nusselt number is found to be in good 
agreement with published correlations. 



Chapter 3 

Data Reduction 


Optical techniques employed in fluid-flow measurements can be classified into 
three categories: (i) direct flow visualization by injecting some tracer particles 
(e.g., dyes, bubbles, solid particles), (ii) laser-Doppler systems where the fre- 
quency shift of scattered illumination from an injected tracer medium in the flow 
field is measured and (iii) measurement of index of refraction of the fluid or its 
spatial derivatives. Three popular methods that fall into the third category are: 

1. Shadowgraph Technique: Here the reduction in light intensity or beam 
divergence is employed. It measures the variation of the second order deriva- 
tive of the index of refraction normal to the light beam. The focus is on 
the reduction of light intensity as the beam diverges. 

2. Schlieren Technique: Here the deflection of light in a variable refractive 
index field is captured. It handles the variation of first order derivative if 
the index of refraction normal to the light beam. Here the focus is on the 
reduction of light in a variable refractive index field is captured. 

3. Interferometric Technique: Here the image formation is related to change 
in the refractive index with respect to a reference environment. It responds 
directly to the refractive index field of the flow system. Here the image 
formation is related to changes in the refractive index with respect to a 
reference environment. 

This chapter discusses the principles of interferometry and procedures for the 
evaluation of the interferograms. The concepts of interferometry with reference 
to temperature field in the infinite fringe setting of Mach-Zehnder interferometer 
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have been discussed. Apart from these features this chapter also discusses the 
techniques used to improve the image quality and fringe thinning. 

3 . 1 Interferometry 

With the infinite fringe setting, the optical path difference between the test and 
the reference beam is zero in the absence of any thermal disturbance. Hence 
interference is constructive and a bright field of view is obtained. In practice, the 
image obtained is fringe-free but may show imperfection associated with the spa- 
tial filter and the interferometer optics. When nonisothermal conditions prevail 
in the path of the test beam (example, a candle flame) each ray of light undergoes 
a change of phase, depending on the extent of change of the refractive index of 
the medium. Hence an optical path difference is established between the test and 
the reference beams, resulting in a fringe pattern. 

As discussed, in interferometry, the image formation is related to the changes 
in the refractive index of the fluid with respect to the reference environment. 
Here attention is restricted to the image patterns that form in the infinite fringe 
setting. When isothermal conditions prevail in the test beam path each ray of 
light undergoes a change of phase, depending on the extent of change , of the 
refractive index of the medium. Hence an optical path difference is established 
between the test and the reference beams, resulting in a fringe pattern. For 
temperature measurements, the refractive index variation has to be related to 
that of temperature. The relationship between the refractive index n and density 
p is 

— — 1 

— 7 — r — = constant (3.1) 

p(n2 -1-2) 

holds for transparent media. For gases, the refractive index is close to unity and 
the expression reduces to 

— — ^ = constant (3.2) 

P 

where p is the density of the fluid. This relationship is called the Gladstone-Dale 
equation (Goldstein, 1983). It follows that dn/dp=constant. Now for moderate 
changes in temperature, density varies linearly with temperature as 


p = p,(l - /5(T - T,)) 


(3.3) 
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Hence dn/dT is also a constant. This implies that changes in temperature 
simultaneously result in changes in refractive index and this leads to changes in 
the phase of the wave. This is the origin of fringe formation in interferometry. 

Let n{r,z) and T{r,z) be the refractive index and temperature field respec- 
tively in the test region. Let rio and To be the reference values of n and T 
respectively as encountered by the reference beam. Let L be the geometrical 
path length covered by the test and reference beams. Now for the test beam, the 
optical length is defined as 


fL 

PL = / n(r, z)ds 
Jo 


(3.4) 


and for the reference beam, 


PL 


=/' 


riods 


(3.5) 


is the optical path length. The interferogram is a fringe pattern arising from the 
optical path difference 


APL= 

Jo 

which in terms of temperature is 

i^PL=^£{T{r,z)-%)ds 


{n{r, z) — no)ds 


(3.6) 


(3.7) 


The integral is evaluated along the path of a light ray given by the coordinate 
s. The fringes seen on the interferograms are locus of points having the same 
optical path difference. Hence on any given fringe the optical path difference 
APL is a constant and 


APL 

(T{r, z) - To)ds = = constant (3.8) 

Hence 

rL 

/ {T{r, z)^s — TqL = constant (3.9) 

Jo 

The integral (T(r,z))ds is defined as TL, where T is the average value of 
T(r, z) over the length L of the laser beam through the test cell. This is also the 
line integral of the function T{r,z). Hence 


L(T — To) = constant 


(3.10) 
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In the infinite fringe setting Equation 3.10 holds good for all the fringes. When 
L is constant for all the rays, Equation 3.10 implies that T is a constant over the 
fringe and hence each fringe represents a locus of points over which the average 
of the temperature field along the direction of the ray is a constant, hence fringes 
are isotherms. 

Defining the function L(T — To) in Equation 3.10 as f{T,L), the fringe tem- 
perature on two successive fringes for same value of L can be given as: 

— APL 

fringe = ^ (3.11) 

fringe 2 :/,(T.L)=(W^ (3.12) 

where A is the wavelength of the laser used. From these two equations, the 
temperature drop per fringe shift, AT^, can be calculated as 

Ar. = i(/,(T.i)-/.(T.i)) = ^, (3.13) 

The number of fringes expected in a projection can be estimated directly from 
the relation 

Number of fringes = (3.14) 

This equality was seen to be satisfied in most of the experiments reported in the 
work. 


3.2 Image Enhancement 

The images of the fringe pattern collected using the CCD camera and the A/D 
card tend to be noisy. Noise can be defined as extraneous and unwanted intensity 
variation which is superimposed over the useful portion of the data. There are 
various factors that are responsible for the generation of noise in the recorded 
data. Examples are the imperfection of the optical components, noise in the 
CCD array and digitization process itself. The most dominant factor is speckle 
which is the major source of noise in the interferometric images. Speckle is a 
super-saturation of the diffraction patterns over the basic interference pattern. 
It arises from the imperfection of the optical surfaces. The presence of micro- 
scopic unevenness on the optical surfaces leads to diffraction. The diffracted rays 





Figure 3.2; Filtered Interferogram 

For Fourier filtering the image is first transformed to the wavenumber space 
with a two dimensional Fourier transform operation. The wavenumber is a kind 


IMAGE ENHANCEMENT 


Figure 3.1; An Original Interferogram 


interfere and generate local fringe patterns and corrupt the global intensity dis- 
tribution. Filtering is required to extract the original intensity information from 
the recorded interferogram which is then processed for contrast improvement. 

The Fourier Filtering method has been used to extract the useful information 
from the recorded interferograms. Figure 3.1 shows the original interferogram as 
recorded from the interferometer. 
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of a spatial frequency and is defined as the number of intensity cycles per pixel. 
Intuitively one can view the low wavenumber harmonics as information and 
high wavenumber harmonics as noise. Now in the wavenumber domain, a two- 
dimensional band pass filter can be used to effectively set the high wavenumber 
components to zero. An inverse Fourier transform restores the filtered image in 
the spatial domain. Figure 3.2 shows the filtered interferogram of the same image 
as shown in Figure 3.1. 

The filtering removes the high wavenumber components of the intensity distri- 
bution and thereby removes the unwanted noise in the image data. At the same 
time, the image becomes blurred owing to the loss of contrast. Since the visual 
appearance of the original interferogram as recorded through the CCD camera is 
not good (as seen from Figure 3.1), the Fourier filtering operation further leads 
to the loss of the contrast. Hence there is a need for contrast improvement of the 
image that has to be analyzed. Different contrast improvement software packages 
are available through which one can enhance the quality of the image. Examples 
are GIMP and Electric Eyes available on the Linux platform. In the present 
work, Electric Eyes software has been used for this purpose. Figure 3.3 shows 
the contrast improved image of the same interferogram as shown in Figures 3.1 
and 3.2. There is a definite improvement in the visual appearance of the image 
and the fringes are now more clearly visible. 



Figure 3.3; Contrast Improved Interferogram 
There are three fringe thinning algorithms suitable for interferograms recorded 
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using the Mach-Zehnder interferometer. In the interferometric images, the fringes 
are represented by the alternate dark and bright bands. The quantitative infor- 
mation (temperature in the present work) is thus spread over the whole width 
of these bands that constitute the fringes. Hence fringe thinning is one of the 
most important operation in the extraction of quantitative information from the 
collected interferograms. It is the process of extraction of the set of points of min- 
imum or maximum intensities in the dark and bright bands of the fringes. The 
minimum intensity corresponds to the points of complete destructive interference 
and hence will have a zero intensity. Similarly, a maximum intensity will appear 
at the points of constructive interference. In the present work, the fringe bands 
are reduced to fringe curves by locating the intensity minima in the dark fringe 
bands. 



Figure 3.4: Superimposed Thinned Image with Original Image 
In the present work, in order to reduce the fringe bands into fringe skeleton 
or curves the following approach has been adopted. This a free-hand drawing 
technique where the midpoints within a dark band of the interferogram are ap- 
proximately located and joined by a smooth curve. It relies exclusively on eye 
judgment. The paint-brush utility of GIMP has been employed in this study. 
The image containing the fringe curves and the original image are subtracted to 
get the fringe skeleton. The superposition of the original interferogram and the 
fringe. skeleton is shown in Figures 3.4 and 3.5. The paint-brush approach has 
the disadvantage of not locating the minimum intensity location, but does not 
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require code development and hence is reasonably fast. 



Figure 3.5: Thinned Image 


3.3 Calculation of Fringe Temperatures 

The calculation of temperature associated with the fringes is the most important 
step in interferometry. Figure 3.6 shows the schematic of a fringe skeleton after 
thinning. The upper and lower walls shown in the figure have known temper- 
atures. Near the walls, a high temperature gradient results in a large number 
of thin fringes which are very closely spaced. Due to the finite resolution of the 
CCD camera, it may be possible that a few number of fringes near the wails could 
be lost during the recording of the interferograms. Hence the first fringe seen on 
the interferogram near the wall will be of some arbitrary order. Therefore one 
cannot assign a temperature to the fringes directly using the value of temperature 
difference per fringe shift AT^ though the wall itself is an isotherm. 

The interferogram is devided into two halves. The upper half is allotted 
temperature on the basis of upper wall temperature and the lower half on the 
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basis of lower wall temperature. 



Figure 3.6: Schematic of an Idealized Fringe Skeleton for Temperature Calculations 
Consider the fringes marked 2,3,4 in Figure 3.6. Fitting a function of type 


T{y) = a + by + cy'^ 

Where y is a vertical coordinate, we get 


(3.15) 


AT, =T 2 -Ts = h{y 2 - yz) + c{yl + yl) (3.16) 

AT, = Ta - T 4 = b{y, - y^) + c{yl + y|) (3.17) 


These two equations are solved for the constants b and c. 

The local wall heat flux temperature is {dT fdy)y-y^ and using Equation 3.15, 
this can be expressed as (6 + 2cyi). Similarly, the gradient in temperature field 
at the first fringe is (6 + 2cy2). Now the average of these two gradients is taken 
and is alloted as the wall temperature gradient. Finally, once this is known, the 
first fringe temperature near the top wall is calculated as 


T, = fe - j,,) + T, 


(3.18) 


The subsequent fringe temperatures are found out by adding or subtracting 
the temperature gradient AT, to the temperature of the previous fringe when 
started from wall. In the present work, these temperature calculations were 
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carried along a particular column. The column was selected where the near wall 
fringes are dense. 


3.4 Nusselt Number Calculations 


Nusselt number calculations have been carried out to determine the heat transfer 
rates at the boundary walls of the cavity, Nusselt number is defined as 


Nu = 


-h dT 


That — Tcold dy 


y=0,h 


(3.19) 


In the present work average Nusselt number at the top and the bottom walls 
have been reported. The average Nusselt number for each of the plates has also 
been compared with the experimental correlation reported by Gebhart et al. 
(1908). For air, the correlation is given by: 


Nu = 1 + 1.44[1 - 


1708, Ra 
Ra ^■^^^5830 


1/3 


- 1 ] 


(3.20) 


For water, the correlation is given by: 


Nu = 1 + 1.44[1 - 1^1 H- - 1] + 2 (^[1 - ^ (3.21) 


The Nusselt number calculations have been performed for all the cases and are 
compared with the reference value as given by the above mentioned correlations. 



Chapter 4 

Results and Discussion 


Experiments performed in the present study are exploratory in nature. They have 
been conducted to examine various issues, such as, coupling between fluid motion 
in air and water, and interfacial effects. To this end, the temperature difference 
across the cavity has been maintained constant at a value of 10°C. Experiments 
have been performed in the following sequence: 

1. Cavity filled with air, 

2. Cavity filled with water, 

3. Cavity half filled with water, with air above. 

All experiments have been continued till steady state is reached. Interfero- 
grams during the transient as well as steady state have been recorded. A steady 
state calculation for the two layer problem that throws light on respective tem- 
perature drops in air and water is presented in Appendix A. 

The present chapter discusses various qualitative and quantitative issues re- 
lated to the evolution of the flow and thermal fields. Flow structures and roll 
pattern, transient evolution of the isotherms and the nature of the steady state 
have been reported. For all the experiments fringe temperature calculations, av- 
erage Nusselt number at the top and the bottom walls and temperature variation 
over the height of the cavity have been presented in the following sections. 

In order to examine the experimental results, a two dimensional numerical sim- 
ulation has been carried out for flow and heat transfer in a square cavity. The pa- 
rameters considered are (i) Ra=34,200, Pr=0.7 and Bi=10.0 and (ii) Ra=10,000, 
Pr=0.05 and Bi=10.0. The governing equations, initial and boundary conditions 
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and the nunaerical technique used along with the results obtained have been dis- 
cussed in Appendix B. 

4.1 Convection in air 

To study convection in air, it is not necessary to employ optical windows, particu- 
larly when the test cell is long in the direction of the light beam. Moreover, quartz 
windows exhibit a refractive index dependence on temperature, that is similar to 
air. Thus difference in the thermal field induced in the optical windows, between 
the test cell and the reference chamber can lead to additional fringes. There is 
also a possibility of a few fringes being lost. In case of liquids, one does not 
have the option of discarding the windows since the liquids has to be confined 
on all sides. The dependence of refractive index of liquids on temperature being 
very weak, in comparison to air as well as quartz, the possibility of distortion 
due to the presence of windows is minimized. In the present section, convection 
experiments in air with and without optical windows have been discussed. 

The experiments in air without optical windows have been considered first. 
The top surface temperature is kept at 16°C and the bottom surface temperature 
at 26'’C. Thus the temperature difference of 10°C is maintained across the cavity. 
This leads to a Rayleigh number of 34,200 in air. The transient evolution of 
the interferograms have been recorded at regular time intervals till the patterns 
become either repetitive or strictly steady in time. 

Figure 4.1(a) and Figure 4.1(b) show the infinite fringe setting of the inter- 
ferometer without and with a candle fiame respectively. The candle flame test 
ensures that wedge fringes are not accidentally present in the field-of-view. Fig- 
ure 4.2 shows the transient evolution of the thermal fleld till the steady state is 
reached. 

The interferograms were collected in regular intervals of time as shown. The 
heat transfer process takes place between the fluid and the surfaces. Since the 
cavity is small, the time required for the temperature of the walls to become 
constant and the flow to reach steady state is also small. For early time, the wall 
temperatures had not stabilized and the corresponding fringe patterns have not 
been shown. With the passage of time, the number of fringes increases till steady 
state is reached. 
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Figure 4.2: 


Evolution Patterns for Air without Windows forTc=16°C, Tft=26°C and 

Ra=34,200 
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The spacing among the fringes near the wall are seen to be small compared to 
the mid-plane of the cavity. The spacing gradually increases towards the fringe 
edges. Since each fringe is an isotherm, a small fringe spacing gives rise to a large 
local heat flux. The largest local heat flux at the cavity walls at steady state 
occurs around the mid-plane of the cavity. The fluid accelerates on one side of 
the mid-plane, reaches a maximum at this point and decelerates to small values 
on the other side of the cavity, as it approaches the side walls. The overall flow 
pattern in the cavity is hence unicellular. 

From Figures 4.2(a-d) it is clear that the isotherm spacings near the top and 
bottom walls of the cavity are less than the spacing in the central region of the 
test section. This is due to the reason that near the boundaries diffusive heat 
transfer rates are higher. In the central region of the test section, the conduction 
heat transfer is practically absent and the dominant mode is convection. 

In the present experiment in a cavity with a square cross-section, the transient 
evolution of the thermal fleld took place over a substantially smaller time scale 
(for about one and a half hour) as compared to cavity square in plan and an 
axisymmetric cavity. Mishra et al. (1998) have reported a time scale of 4-5 hours 
in their experiments. Once the flow fleld has fully evolved, the fringe patterns 
were observed to be quite steady without any sideways movements of the fringes. 
This sideways movement of fringes was observed for higher Rayleigh numbers. 
Figure 4.2(b), Figure 4.2(c) and Figure 4.2(d) show the truly steady state nature 
of the flow fleld. These interferograms were recorded after one and half hours of 
the experimental run time to ensure the steady state behavior of the convection 
pattern. The number of fringes in the present experiment are more than that 
estimated by the formula N = {Th — T^jhT^. Some false fringes are captured in 
the central region of the cavity where the temperature gradients are small. 

Figure 4.3(a) and Figure 4.3(b) show the original steady state image and 
the Fourier Altered image respectively. Figure 4.3(c) shows thinned image and 
Figure 4.3(d) shows superimposed thinned image. From the thinned image, the 
calculations for the temperature of the fringes and the local Nusselt number can 
be performed. Figure 4.3(c) shows that the fringe thinning algorithm has been 
successful in the present work. 
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Figure 4.3: (a) Original Interferogram (b) Filtered Interferogram 



(c) (d). 

Figure 4.3: (c) Thinned Image (d) Superimposed Thinned Image 
for Air when the cavity is filled with Air for Ra=34,200 
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(a) (b) 

Figure 4.4; Temperature Profiles in air filled cavity for Ra=34,200 
at (a) x=1.07 mm and (b) x=21.4 mm 



X 


Figure 4.5: Local Nusselt number variation at the top and bottom walls 
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At steady state, the variation of the temperature field is plotted as a function 
of the vertical distance as shown in Figure 4.4 for two equally spaced columns 
along the width of the cavity. 

Heat transfer rates at the boundary walls in terms of the local Nusselt number 
can be calculated by using Equation 3.19 as 


Nu = 


-h dT 


That — Tcoid dy 


y—0,h 


Figure 4.5 shows the variation of local Nusselt number at the top and bottom 
walls of the cavity. ‘ The individual plate-averaged Nusselt numbers for the case 
of air without windows are 3.18 and 3.12 at the cold and the hot surfaces respec- 
tively of the cavity. The average Nusselt number for each of the plates has also 
been compared with the experimental correlation given by Gebhart et al. (1988), 
Equation 3.20. The value of average Nusselt number as calculated from Equa- 
tion 3.20 is 3.17 for this case. The comparison between the present experiments 
and the published correlation is thus seen to be quite good. 

Analysis of the flow field observed in air in the presence of optical windows 
is now presented. The experimental conditions are identical to the previous case, 
namely, Tc=16°C, T/i=26°C, AT=10°C and Ra=34,200. The cavity is enclosed 
by optical windows. Air is thus confined on all sides and contact with the out- 
side atmosphere is prevented. The transient evolution of the flow field has been 
recorded at regular time intervals till the steady state is reached. Here the num- 
ber of fringes is lower than in the experiments without optical windows. However 
the number of fringes match the estimate based on AT^. The fringe patterns 
are qualitatively similar to those without windows. The number of fringes were 
found to be 7 at steady state. 
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Figure 4.6: Infinite Fringe setting (a) without and (b) with Candle Flame 
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Figure 4.7: Evolution Patterns for Air with Optical Windows for Tc=16°C, T/i=26°C and 

Ra=.^4 9nn 
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Figure 4.6 shows the infinite fringe setting without and with a candle flame in 
the presence of optical windows. Figure 4.7 shows the transient evolution of the 
thermal field till the steady state is reached. The interferograms were collected 
in regular intervals of time as shown. The presence of windows did not alter the 
time required to reach steady state. With the passage of time, the number of 
fringes increased till the steady state was reached. 

Once the flow field was fully evolved, the fringe patterns were observed to be 
quite steady and motionless. Figure 4.7(b), Figure 4.7(c) and Figure 4.7(d) show 
the steady state nature of the flow field. These interferograms were recorded 
after one and half hours of the experimental run time to ensure the steady state 
behavior of the convection pattern. 

Figure 4.8(a) and Figure 4.8(b) show the original steady state image along 
with the Fourier filtered image respectively. Figure 4.8(c) shows thinned image 
and Figure 4.8(d) shows superimposed thinned image for air in the presence of 
optical windows. From the thinned image, calculation for temperature of fringes 
and the local Nusselt number can be carried out. 

Figure 4.9(a) and Figure 4.9(b) are the Fourier filtered images of Figure 4.7 
(d) with different filtering parameters and Figure 4.9(c) and Figure 4.9(d) are 
the thinned image and the superimposed thinned image respectively. From the 
calculation of temperature of fringes and Nusselt number, it was observed that 
there was no effect of the filtering parameter on data analysis. 

At steady state, the variation of the temperature field plotted as a function 
of the vertical distance is shown in Figure 4.10 for two equally spaced columns 
of the cavity. It is to be recalled that the y-coordinate is measured from the cold 
top wall. 

Heat transfer rates at the boundary walls in terms of the local Nusselt number 
has been calculated by using the Equation 3.19. 




Figure 4.8: (c) Thinned Image (d) Superimposed Thinned Image for air when 
the cavity is with optical windows for A.r for Ra=34,200 
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Figure 4.10: Temperature Profiles for two equally spaced columns 
of cavity for (a) x=1.07 mm and (b) x=21.4 mm 
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Figure 4.11: Local Nusselt number variation at the top and battom walls 
for Air when the cavity is with optical windows for Ra=34,200 
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Figure 4.11 shows the variation of local Nusselt number at the top and bottom 
walls of the cavity. The individual plate averaged Nusselt number for the case of 
air with windows were found to be 2.90 and 3.23 at the cold and the hot surfaces 
respectively of the cavity. The average Nusselt number for each of the plates has 
also been compared with the experimental correlation given by Gebhart et al. 
(1988), Equation 3.20. The value of average Nusselt number as calculated from 
Equation 3.20 is 3.17. The Nusselt number calculated in the present experiment 
matches well with the reference correlation. 

4.2 Convection in a water filled cavity 

Experiments with the test cell filled with water are reported in the present sec- 
tion. The hot and the cold walls of the cavity were maintained at 16 and 26‘’C 
respectively, the effective temperature difference is 10°C. to get meaningful fringe 
patterns, it was found necessary to fill the reference chamber also with water. 
The Rayleigh number for the present experiment was calculated to be 16,69,314. 
Interferograms were recorded at various time instants, but no clear steady state 
was to be seen. This is understandable because the Rayleigh number for the 
present experiment stretches well into the turbulent regime. Keeping in mind the 
low sensitivity of refractive index of water to temperature, one can anticipate a 
large number of fringes in water. 

Figure 4.12 shows the infinite fringe setting of the interferometer without and 
with a candle fiame. The fringes for the candle flame in the presence of the test 
cell with water confirm that the optical alignment has been proper. Figure 4.13 
shows the transient evolution of the thermal field. Since the Rayleigh number 
is extremely high, the fringe patterns are seen to change rapidly. The patterns 
show that the flow is in three-dimensional turbulence zone. 

Figure 4.14(a) and Figure 4.14(b) show the original image along with Fourier 
filtered image respectively. 

Figure 4.14(c) shows thinned image and Figure 4.14(d) shows superimposed 
thinned image for a cold wall temperature Tc=16 "C and a hot wall temperature 
Tft= 26 °C. 
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Figure 4.14: (c) Thinned Image (d) Superimposed Thinned Image for Water Ra— 16,69,314 
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In the present experiment it was not possible to clearly identify the near wall 
fringes for this reason the local as well as the average wall heat transfer rate 
could not be calculated. Continuous fringes however could be seen in the middle 
of the test cell. Transverse heat transfer based on the fringe spacing near the core 
showed the Nusselt number to be close to 0.07. This is to be contrasted with the 
correlation of Gebhart et al. (1988) Equation 3.19, which gives a Nusselt number 
of 10.0. The difference between the two arises from the fact that the wall heat 
transfer is diffusion dominated, while in the core of the cavity heat transfer is 
predominantly by advection. 

4.3 Convection in a cavity half filled with water 

Experiments have been conducted with the test cell filled with water upto the 
raid-plane. The cavity thus contains two fluid phases, namely air and water. 
To generate fringe patterns successfully with the interferometer, the reference 
chamber was also maintained with water upto its mid-plane. An unexpected 
difficulty was encountered during the alignment of the test and reference beams 
with respect to the interface. Owing to mild refraction eflfects, the interface 
appeared on the screen as a dark band. The width of the band was brought to 
within limits by improving the alignment of the apparatus repeatedly. Thermal 
fields above and below the interface could be recorded by the interferometer and 
being seen to be insensitive to the width of the band. It was however found 
necessary to check for the alignment of the interferometer individually for the air 
and the water phases. 

Experiments have been conducted with the hot and cold surfaces maintained 
at 26 and 16"C. The corresponding Rayleigh number in air and water were cal- 
culated to be 4,185 and 16,413 respectively. These numbers are based on the half 
cavity height and the interface temperature calculated as in Appendix A. 

4.3.1 Convection in the air phase 

Figure 4.15(a) and Figure 4.15(b) show the infinite fringe setting of the interfer- 
ometer without and with the candle flame respectively. The black band in the 
center is due to the interference of the reflected and the refracted beams of light 
from the water and air portions of the cavity. The fluid above the band is air and 
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Figure 4.15: Infinite Fringe setting (a) without and (b) with Candle Flame 
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c) 2:00 

4.16: Evolution Patterns for air when the cavity is half filled with Water for Tc=16‘=C 
Anri Ra ('a.ir'i=4.185 and Ra(water)=16,4l3 










convection in A CAVITY HALF FILLED WITH WATER 


Figure 4.17: (c) Thinned Image (d) Superimposed Thinned Image for Air for Tc-16°C 
T/i=26°C and Raa^r— 4,185 and Rau;ater— 15)413 
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below is water. 

Figure 4.16 shows the transient evolution of the thermal field in air till the 
steady state is reached. The interferograms were collected at regular intervals of 
time. In this case also steady state obtained very fast. The fringes are not straight 
lines but have some curvature. From the figure it is clear that the isotherm 
spacings near the top wall and the interface are less than the spacing in the 
central region. Here the number of fringes is lower than that of the experiment in 
which the cavity is filled with air. This is because the Rayleigh number is lower. 
In the present experiment, we encountered a new problem of dew formation on 
the optical windows. Hence ail the fringes could not be captured. Reflection of 
light from water in the air portion of the cavity is seen in Figure 4.16. 

Once the flow fleld was fully evolved, the fringe patterns were observed to be 
quite steady. Figure 4.16 shows the steady state nature of the flow field. These 
interferograms were recorded after one and half hours of the experimental run 
time to ensure the steady state behavior of the convection pattern. The number 
of fringes is around 3 at steady state. 

Figure 4.17(a) and Figure 4.17(b) show the original steady state image along 
with Fourier filtered image respectively. Figure 4.17(c) shows thinned image 
and Figure 4.17(d) superimposed thinned image for air for Ra=4185. From the 
thinned image, local temperature and Nusselt number are calculated. 

4.3.2 Convection in the water phase 

For getting the fringe patterns in water, the water layer in the test cavity was 
matched with that of the reference cell. Infinite fringe setting was then achieved 
for the water phase alone. 

Figure 4.18(a) and Figure 4.18(b) show the infinite fringe setting of the in- 
terferometer without and with the candle flame respectively. The black patch 
once again represents the interface. Figure 4.19 shows the transient evolution of 
the thermal field in water till the steady state is reached. The interferograms 
were collected in regular intervals of time. Figure 4.19(a) has fewer number of 
fringes than at steady state, namely Figure 4.19(c) and Figure 4.19(d). Figure 
4.19(b) was recorded just before the steady state reached. The number of fringes 
increases continuously till steady state reached. 






Interface 


Infinite Fringe setting (a) without and (b) with Candle Flame 


Figure 4.19: Evolution Patterns for Water when the cavity is half filled with Water for Tc-16°C 
Th=26°C and Raair=4,185 and Ra„,ater=16,413 
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Figure 4.18: 
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Hot Wall 


Figure 4.21; 




(a) Raw image with extraneous features and (b) Modified image 
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Interface 


Figure 4.23: Temperature Profiles in the cavity half filled 
with water at (a) x=1.07 mm and (b) x=21.4 mm 
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(a) Mechanical Coupling 

Figure 4.24: (a) Local Nusselt number variation at the top and 
bottom walls (b) Rolls formation in air and water 
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The fringes are not straight lines but have a curvature. From the figure it is 
clear that the isotherm spacings near the interface and bottom wall of the cavity 
IS less than the spacing in the central region of the test section. Here the number 
of fringes is lower than for the experiment in which the cavity is fully filled with 
water. This is because the Rayleigh number is lower. Once the flow field has fully 
evolved, the fringe patterns were observed to be quite steady. Figure 4.19 shows 
the steady state nature of the flow field. These interferograms were recorded 
afi.ei one and half hours of the experimental run time to ensure the steady state 
behavior of the convection pattern. The number of fringes is around 12 at steady 
state. 


Figures 4.20(a) and Figure 4.20(b) shows the original steady state image along 
with Fourier filtered image respectively. Figure 4.20(c) shows thinned image and 
Figure 4.20(d) shows the superimposed thinned image at Ra=16,413. From this 
thiniKxl image!, calculation for temperature of the fringes and Nusselt number can 
be carried out.. 

While working with two fluid layers of air and water, unexpected difficulties 
were experienced during interferometry. Two such factors already referred above 
are the liroadening of the interface into a black patch and dew formation on optical 
windows. In addition, mild refraction effects from the cavity walls got amplified 
over the balancing arms of the interferometer. The image thus formed was larger 
than the true size of the cavity. Certain ghost features were also revealed in the 
interferogram. It was however possible to eliminate the undesirable features of the 
image by using an eye-judgement. The original and the modified images showing 
the extraneous features as well as their removal are shown in Figures 4.21(a) 
and (1)) respectively. The Figure 4.22(a) and Figure 4.22(b) are the steady state 
pattern for air and water when the cavity is half filled with water and air. 


At steady state, the variation of the temperature field plotted as a function 
of the vertical distance as shown in Figure 4.23 for two equally spaced columns 
along the width of the cavity. The y-coordinate is measured from the interface. 
From the figure it is seen that there is more temperature drop in air than in 
water. Figure 4.24(a) shows the local Nusselt number variation at the top and 

bottom walls of the cavity. 

Figure 4.24(b) shows the roll formation in both the portions. From the convec- 


tion pattern it is understai^ab.le 


rs are mechanically coupled. 
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Heat transfer rates at the hot and the cold walls in terms of the local Nusselt 
number have been calculated by using Equation 3.19. Figure 4.24(a) shows the 
variation of local Nusselt number at the top and bottom walls of the cavity. The 
individual plate averaged Nusselt number for the case of two layer problem are 
1.26 and 2.33 at the cold and the hot surfaces respectively. The average Nusselt 
number for each of the plates has also been compared with the experimental cor- 
relation given by Gebhart et al. (1988), Equations 3.20 and 3.21. The value of 
aveiage Nusselt number as calculated from Equation 3.20 is 1.85 for air and from 
Equation 3.21 is 3.68 for water. There is a difference between the Nusselt number 
from the present experiment and the correlation. The correlation considers the 
temperature of hot and cold walls where the heat transfer mechanism is conduc- 
tion. In the present experiment one of the boundaries is an interface between the 
two fhiids and not rigid wall. Heat transfer at the fluid interface is dominated by 
conv(Hd,ive effects. 

4.3.3 Experiments with a new set of surface temperatures 

For most fluids, all thermophysical properties including thermal conductivity and 
viscosity are temperature dependent. To examine the sensitivity of the thermal 
field to the average temperature level in the cavity, experiments were repeated 
with new temperature levels of 34 and 24‘’C, keeping the cavity temperature 
difference at the original value of 10°C. The original fringes in air and water as 
well as the processed images are presented in Figure (4.25-4.31). These figures 
are seen to be remarkably similar to those presented in Sections 4.3.1 and 4.3.2. 

The (essential conclusion to emerge from the above experiment is that the ther- 
mal field is buoyancy dominated, being dependent on the changes in density with 
the changes in temperature. The temperature dependence of the conductivity 
and viscosity play a very insignificant role in deciding the convection pattern. 
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Water 


(a) (b) 

Figure 4.25: Infinite Fringe setting (a) without and (b) with Candle Flame 


1*45 (c) 

Figure 4.26: Evolution Patterns for Air when the cavity is half filled with Water for Tc-24 C, 
T,,=34‘’C, Raair=3,699 and Ra^ater=22,350 
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Figure 4.29: Evolution Patterns for Water when the cavity is half hlled wih Water tor T 

Th=U°C, Ra«r=3,699 and Ra^ater=22,350 
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Figure 4.28: Infinite Fringe setting (a) without and (b) with Candle Flame 


24°C, 








A CAVITY HALF FILLED WITH WATER 





Interface 



' — — I Water 

Thermal Coupling 

Figure 4.31: R,oll patterns in air and water phase for Raa, >=3,699 and Ra^ater=22,350 

A f(iw minor differences were however recorded with the experiments employ- 
ing higlKir surface temperature levels. These are: 

1. The fringe pattern in the air phase shows a reversal in the direction of 
circulation, while that in water has remained unchanged. This obstruction 
weakens the result of mechanical coupling between the two phases clearly 
recorded in the experiments at a lower temperature. 

2. The air water interface is not as clearly delineated when the average cavity 
temperature is higher {29°C, as in the experiments with a lower average 
cavity temperature of 21°C). This can be explained partially by the reduc- 
tion in surface temperature with increase in temperature and partially by 
the absence of a shared interface between the test and the reference beam. 



Chapter 5 


Coiiclnsioiis and Scope for Future 
Work 


convection in two-dimensional differentially heated fluid layers 
ha,s IxM'n (‘X|)(*rini('ntally stiidied using Mach-Zehnder interferometer. The cavity 
ha,s InH'ii ta.k('n to he (a) filled with air, (b) filled with water and (c) half filled 
wifii wa.t,(vr. 'I’Ik' driving potential in all the experiments was 10°C. The following 
conclusions W('r(^ arrived at in the present study: 

1. A st.('a,(ly two-dimensional convection pattern in air was seen to be estab- 
lished. rh(^ average Nusselt number for the cavity matched well with the 
pul)lished correlation. 

2. Wlien tlu^ cavity is filled with water, the Rayleigh number was extremely 
liigh and the flow was unsteady and turbulent. 

3. wIk'h th(' cavity was half filled with water, the interface temperature was 
sueli a.s to produce a lower Rayleigh number in both air and water. Nearly 
st,('a.dy convoH'.tive field were seen to establish in both the phases. The 
int,('rra,c(^ t,emperature calculated from the water side was found to be quite 
clos(! (,o that from the air side. For a cavity average temperature of 21"C, 
th(' t,wo laycirs were seen to be mechanically coupled. For a cavity average 
l,(‘nii)(!rature of 29"C, the two layers were thermally coupled. 

'I'lu^ i)r(‘s<'nt research can be extended in the following manner to cover other 
areas of pr'act;ical interest. 

1. Nffeci, of layer height on the nature of coupling between air and water. 


- ^ 

2. Convection in low Prandtl number liquid with an overlaid air cushion. 

3. Convection in superimposed immiscible liquid layers. 
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Appendix A 

Numerical Simulation for 
Interface Temperature 


A.l Design of Test Cell 


Tins apixuitlix throws light on the steady state calculation for the respective 
tenii){'ratur(^ drops in air and water for the two layer problem. 

h'or t,h(' c:onv(H;t,iv(! How to be in 2-dimensional steady state regime for air and 
water in a t,wo plia.s(> convection problem the Rayleigh numbers should fall in the 
flow iH^gimt! given by if.. Krishnamurti (1970), Figure 2.2. Individual Rayleigh 
numbers ar<^ calcuhited by using equation 


av 


Ra = 


g^ATd^p\ 

flK 


(A.1) 

(A.2) 


wh(‘r (' d is th(‘ half of the cavity height. 

h'or this ll(u-(g th(‘ Nusselt numbers for air and water are calculated by using 
correflations given by Hollands et al (1975). For air, the correlation is given by: 


Nua = 1 + 1.44[1 


Ra J ^^5830^ ^ 


(A.3) 


And for water, the correlation is 
Nil,,, 1 ) 1.4'l[l 


+ + (A.4) 


140 


140 


Ra ^ ' ^^5830 

Where the terms within square brackets are taken as zero if they have negative 
value. The additional term for water is suggested by Hollands et al. (1975) to be 



A.l. DKSKJN OF TEST CELL 



(hie t,o iM ('xistencx^ of thermals. For Ra< 1708 
oc,cvirs by eoiiduetiou and Nu=l. 


the heat across the fluid layers 


By IJw law i>l (aniservation of energy in terms of heat given by 
heat <am<l„<a«i through bottom plate=heat convected in water 
=hcat cOTivcctcd in air=heat conducted in top plate 

It is assumed that the thicknese of the top and bottom surfaces are much smaller 

than that of the finid layers. Hence temperature gradients occur only due to 

convcKd.ion in fluids. Using 


wlmv 


and 


lu'at convected in water=heat convected in air 

hn,{Th - Ti)A = ha{Ti ~ Tc)A 

NUio Rft/r 


Jhn 


hn 


NUaRfa 


(A.5) 

(A.6) 

(A.7) 


whm'c (i,„ and d„ a.r(' the charecteristic length of the fluids. In the case of cavity 
filled witch half wabvr charecteristic lengths are equal = da)- Hence 


Nu,„A'„,(T„ - Ti) = l^naKaiTi - %) 


(A.8) 


Th(i inti'rfaei' t.einjxu’ature Tj is calculated from the above equation. For a given 
cavity Inuglit a.nd top and bottom wall temperatures, one can get the properties 
of th(‘ fluids, since the properties of the fluids such as p, Cp, p, K, a, v, depend 
on t,(‘mp('ra.l.ur('. Curv(^ fitting is done by using equation: 


y = do A diT + 02 ^^^ 


(A.9) 


when* ■// can fa* p, Cp, p, K, jd, a, and v. For air ^~Taverage’ 
propi'rtif'K of air and water are obtained for the required temperature. For air 
th(^ prop(‘rl,ies ar(^ taken at the average temperature For water the average 

temperature is wluire is the cold (upper) wall temperature and Th is the 

hot (low'er) wall teuqxu'ature, Ti is the interface temperature. 

Al’tf'r ea.lculating the interface temperature, number of fringes in water and 

air can be ealeulai.ed by using formula 


muuber of fringes in water = 


{Th-Tj) 

^^£(water) 


(A.10) 
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and 


mmiber of fringes in air = 


{Ti-n) 


AT 


£{air) 


(ATI) 


resixH^l.ivt'ly. 

Using a,l)c)V(‘ formulation for different cavity heights and for different surface 
tenip('i<itm('s llayleigh numbers, Prandtl numbers, interface temperatures are 
cal(mla.t('d for water and air, which are shown in Table-A.l, Table-A.2 and Table- 


A.3. 


Table-A.l 


h 

I'll 

Tc 

Ra (air) 

Ra (water) 

Nu(hot) 

Nu(cold) 

Ti 


•li\ 

10 

1116 

2538 

2.120 

1.0 

19.81 

20 

M) 

20 ! 

954 

3907 

2.527 

1.0 

29.84 

mm 

35 

15 

1913 

8992 

3.186 

1.155 

34.70 


Table-A.2 


h 

I'll 

Vc 

Ra (air) 

Ra (water) 

Nu (hot) 

Nu (cold) 

Ti 


l'" 2() 

10 

2179 

5085 

2.717 

1.311 

19.80 

25 

.30 

20 

1865 

7215 

2.996 

1.121 

29.85 

mm 

' ' 35 

15 

3725 

21830 

3.923 

1.779 ^ 

34.63 


Table-A.3 


h 

“Th''^ 1 

Tc " ‘ 

Ra (air) 

Ra (water) 

Nu (hot) 

Nu (cold) 

Ti 


' 20 

10 

3756 

9983 

3.265 

1.785 

19.78 

30 

30 

20 

3212 

15350^ 

3.629 1 

1.674 n 

29.81 

mm 

■' '35 

’ 25 

2975 

18177 

3.769 

1.614 

34.82 


b'roin t.he tabhw it is observed that the temperature drop in air is much higher 


than tliat of water. 



Appendix B 

Numerical Simulation of Fluid 
Flow and Heat Transfer 


1 h(^ pr(!S(^nt. appendix is concerned with the numerical solution of the buoyancy- 
driven flow in a fluid fllled square cavity. The dimensions of the test cell are 
tilu^ sanui as considered in the experimental work. The numerical solutions have 
hi'Mi i)resented for Ra=34,200 and Pr=0.7 (air) and for Ra=10,000 and Pr=0.05 
(liquid nHd,al). Open as well as closed cavities have been considered. 

Natural convection in the fluid is governed by the laws of conservation of mass, 
momentum and energy. These are mathematically represented by the continuity 
equation, Navier-Stokes equation and the energy equation, which are represented 
in 2-dim(msional cartesian system as follows: 


Continuity: 

;r-moui(MitiUm: 

;(/-mom<ntum: 

IDnergy: 


du _ Q 

dx dy 

. du du dp r 


u 


p[u^ + 

^ dx dy' dy 

. dT , dT, 


(B.l) 

(B.2) 

(B.3) 

(B.4) 


Th(^ a,l)ov(^ (filiations have been recast in a nondimensional form using the follow- 
ing quantities: 


Dimensionless temperature: 


0 = 


T-Tc 

Th-Tc 




(B,5) 
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Dimensionless length: 


Dimensionless velocity: 


(- y 

h 

(u,v 


vjh 


(B.6) 

(B.7) 


Dimensionless pressure: 

p = y p°yy 

p{p/hy 


(B.8) 


The non-dimensional parameter responsible for the buoyancy induced flows in 
the cavity namely the GrashofF number is defined as 


Gr = ^ 
Pr 


(B.9) 


where R,a is the Rayleigh number and is defined as 


ua 

(B.IO) 

and Pr is the Prandtl number defined as 


Pr = - 
a 

(B.ll) 

Biot number Bi is defined as 

k 

(B.12) 

Marangoni number Ma is defined as 


Ma= 

pa 

(B.13) 

To solve the governing equations for flow and temperature, a 'ip— co- 
hos been employed. The u and v velocities have been defined as 

-6 formulation 

dtp 

(B.14) 

dtp 

^ dx 

(B.15) 

One can now derive the following: Stream Function Equation: 


dv du 

Since, 0 ; = ^-;^ 

(B.16) 

d'^tp d'^tp 

(B.17) 
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Vorticity Equation: 


Energy Equation: 


du) duj ,Ra.50 

u \- V — = ( 1 1 4 - ( i 

dx dy ^dx'^ dy'^’ ^Fx’dx 


de 86 1 ^d'^e 8^, 

“3r + ”iir = 5:te + ^) 


(B.18) 


(B.19) 


dx dy Pr'dx^ dy^' ' ’ 

For computation of the stream function, vorticity and temperature, three 
geometries have been studied: (1) cavity is closed, (2) open from the top, and (3) 
the surface tension term is included at the free surface. The boundary conditions 
for the three sides i.e., for the lower wall and the two side walls are identical for 
all the three cases. But for the top walls the boundary conditions are different. 
The top wall boimdary conditions for all the cases can be represented as 
CLOSED CAVITY: 

Ip = 0 (B.20) 

.. = (B.21) 

8y^ 

-~=Bie (B.22) 

dy 


(B.22) 


OPEN CAVITY: 

i/> = 0 

w = 0 
dy 

OPEN CAVITY with SURFACE TENSION: 


(B.23) 

(B.24) 

(B.25) 



xP = 0 

(B.26) 


A. 

00 = Ma— 

8x 

(B.27) 


-^=BM 

dy 

(B.28) 

The remaining 

boundary conditions are as follows: 


Bottom Wall: 

i/) = 0 

(B.29) 


8‘^'ip 

(B.30) 
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e = i 

(B.31) 

Side Walls: 


-ij; = Q 

(B.32) 

u = 

dx^ 

(B.33) 

de 

dx~^ 

(B.34) 

At All Corners: 


d^'ih 

2-0 ^^2 

(B.35) 

The following initial conditions have been applied: 


Ip = 0 

(B.36) 

w = 0 

(B.37) 

0 =l-y 

(B.38) 

The governing equations have been discretised 

using an implicit finite differ- 

ence scheme. A uniform grid of 81x81 is used. 

The convective terms of the 


momentum equations are discretized using the 2nd order upwind scheme while 
central difference discretisation has been applied to the diffusive terms. The ma- 
trix inversion is by the Gauss-Siedel method. The results of numerical simulation 
for Ra=34,200 and Pr=0.7 and for Ra=10,000 and Pr=0.05 have been presented 
below. For these two fluids three cases are considered: 


1. Cavity is closed, 

2. Cavity is open at the top surface, 

3. Cavity is open, but surface-tension on the free surface is non zero. 

Figure B. 1(a), Figure B.l(b), Figure B.l(c) and Figure B.l(d) show the nu- 
merically generated velocity vectors, streamlines, isovorticities and isotherms for 
Ra=34,200 and Pr=0.7 for a closed cavity. As seen from the plots of the isotherms, 
the fluid accelerates on one side of the mid-plane, reaches a maximum at this point 
and decelerates to small values on the other side of the cavity, as it approaches 
the side walls. The overall flow pattern in the cavity is hence unicellular with the 
roll moving in a clockwise direction. 
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Figure B.4: Temperature Profiles at 3 columns in (a) Closed cavity 
(b) Open cavity without and (c) with Surface-Tension at x=0.25, 
0.50, 0.75 for Ra=34,200, Pr=0.70, Bi=10.0 and Ma=1.0 









B.l. NUMERICAL SIMULATION FOR RA=10,000 AND PR=0.05 
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Figure B. 2(a), Figure B.2(b), Figure B.2(c) and Figure B.2(d) show the nu- 
merically generated velocity vectors, streamlines, isovorticities and isotherms for 
Ra— 34,200 and Pr=0.7 for an open cavity without considering the effect of sur- 
face tension. The overall flow pattern in the cavity is unicellular with the roll 
moving in a clockwise direction. 

A value of Pr=0.7 is generally seen for a fluid medium in the gaseous phase. 
If the working fluid is a gas, free surface boundary conditions as discussed in 
the present Appendix are not applicable. There is partial evidence that the flow 
patterns in gases as well as liquids are quite similar, particularly at lower Rayleigh 
numbers. For this reason, the flow fields computed with Pr=0.7 can be expected 
to be structurally similar to those in a liquid such as water, whose Prandtl number 
is in the range of 5-10. The above analogy can not be extended to low Prandtl 
number fluids such as liquid metals. 

Figure B. 3(a), Figure B.3(b), Figure B.3(c) and Figure B.3(d) show the nu- 
merically generated velocity vectors, streamlines, isovorticities and isotherms for 
Ra=34,200 and Pr=0.7 for an open cavity including the effect of surface tension. 
The overall flow pattern in the cavity is unicellular with the roll moving in a 
clockwise direction. Figure B.4 shows the temperature variation as a function of 
vertical coordinates at three different equally spaced columns along the length of 
the cavity for all the three cases disscussed above. 

B.l Numerical Simulation for Ra=10,000 and 
Pr=0.05 

The results of numerical simulation for Ra=10,000, Pr=0.05 and Bi=10.0 have 
been presented in this section. It may be noted that a low prandtl number is 
encountered with liquid metals. 

Figure B.5(a), Figure B.5(b), Figure B.5(c) and Figure B.5(d) show the nu- 
merically generated velocity vectors, streamlines, isovorticities and isotherms for 
Ra=10,000, Pr=0.05 and Bi=10.0 for a closed cavity. As seen from the plots of 
the isotherms, the fluid near the walls accelerates on one side of the mid-plane, 
reaches a maximum at this point and decelerates to small values on the other 
side of the cavity, as it approaches the side walls. The overall flow pattern in the 
cavity is hence unicellular with the roll moving in an anticlockwise direction. 











B.l. NUMERICAL SIMULATION FOR RA=10,000 AND PR=0.05 


Figure B.7: Numerically generated (a) Velocity Vectors, (b) Streamlines, (c) Isovorticities 
and (d) Isotherms for Ra=10,000 and Pr=0.05 for an Open cavity with Surface-Tension 



B.l. NUMERICAL SIMULATION FOR RA= 10,000 AND PR=0.05 
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Figure B.8: Temperature Profiles at 3 columns in (a) Closed cavity 
(b) Open cavity without and (c) with Surface-Tension at x=0.25, 
0.50, 0.75 for Ra=10,000, Pr=0.05, Bi=10.0 and Ma=1.0 






B 2. effect of SURFACE-TENSION ON THE SHAPE OF THE FREE 
SURFACE 


However smaller cells are also visible here near the top right and lower left 
corners. 

Figure B. 6(a), Figure B.6(b), Figure B.6(c) and Figure B.6(d) show' the nu- 
merically generated velocity vectors, streamlines, isovorticities and isotherms for 
Fta=10,000 and Pr=0.05, Bi=10.0 for an open cavity without considering the ef- 
fect of surface tension. The overall flow pattern in the cavity is bicellular with 
two rolls moving in clockwise and anticlockwise directions. 

Figure B. 7(a), Figure B.7(b), Figure B.7(c) and Figure B.7(d) show the nu- 
merically generated velocity vectors, streamlines, isovorticities and isotherms for 
Ra=10,000, Pr=0.05 and Bi=10.0 for an open cavity including the effect of sur- 
face tension. The overall flow pattern in the cavity is bicellular with two rolls 
moving in clockwise and anticlockwise directions. The strength and the size of 
the rolls have now changed in relation to Figures B.6(a-d), the low'er roll being 
more prominent. 

Figure B.8 shows the temperature profiles for liquid metal for closed, open 
without and with surface-tension for equally spaced three columns. 

B.2 Effect of Surface-Tension on the shape of 
the Free Surface 

Figure B.9(a) and Figure B.9(b) show the surface deformation due to surface- 
tension for a fluid with Pr=7.0, a cavity height /i=16.05mm with ^=0, that is 
for constant a and for ^=constant that is for nonconstant a respectively. For 
water, the contact angle at the cavity walls was taken to be 0.78 radians. The 
figures show that the displacement of free surface owing to convection is 5.4 mm 
for constant surface-tension and 6.9 mm for non-constant surface-tension. Figure 
B. 10(a) and Figure B. 10(b) show the surface deformation due to surface-tension 
for a fluid with Pr=0.05, a cavity height /i=16.05 mm with ^=0, that is for 
constant a and for ^=constant that is for nonconstant a respectively. For a 
liquid metal, the contact angle at the cavity wall was taken to be 2.34 radians. 
The calculations show that the displacement of firee surface owing to convection is 
6.9 mm for constant surface-tension and 7.0 mm for non-constant surface-tension. 



B.2. EFFECT OF SURFACE-TENSION ON THE SHAPE OF THE FREE 
SURFACE 
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Figure B.9(a): Effect of Constant Surface-Tension on the displacement of the free 
surface owing to convection in water, Ra=34,200, Pr=7.0, Bi=10.0, Ma =1.0 
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Figure B.9(b): Effect of Non-Constant Surface-Tension on the displacement of the free 
surface owing to convection in water, Ra=34,200, Pr=7.0, Bi=10.0, Ma =1.0 
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Figure B.lO(a): Effect of Constant Surface-Tension on the displacemen t of the free 
surface owing to convection in a liwuid metal, Ra=10,000, Pr=0.05, B i=10.0, Ma=1.0 
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Figure B. 10(b); Effect of Non-Constant Surface-Tension on the displacement of the free 
surface owing to convection in a liquid metal, Ra=10,000, Pr=0.05, Bi=10.0, Ma=1.0 



